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Z: [H5] EZ] 2-78 83 PESCE, Culex 种 ( 亚 种 ) 的 CO 工序 列 和 ITS2 序列 进行 测序 ,构建 


和 讨论 这 些 种 ( 亚 种 ) 的 分 子 系统 发 育 关系 。【 方 法】 测定 了 库 蚊 属 20 个 种 ( 亚 种 ) 的 CO T de ITS2 
序列 ,并 从 NCBI 数据 库 中 下 载 了 库 蚊 属 另外 20 个 种 ( 亚 种 ) 的 COT 和 3 个 种 ( 亚 种 ) 的 ITS2 序 
列 。 对 库 蚊 属 40 个 种 ( 亚 种 ) 的 CO 工序 列 和 其 中 23 个 种 ( 亚 种 ) 的 ITS2 序列 进行 碱 基 构 成 .种 间 
遗传 距离 和 饱和 度 分 析 , 并 对 CO D. + ITS2 序列 进行 ILD (incongruence length difference) 检验 。 分 
别 使 用 2 种 分 子 数据 集 (CO D 和 ITS2) 的 核 普 酸 序列 ,用 最 大 似 然 法 ( ML)、 贝 叶 斯 法 ( BI) 、 邻 接 法 
(NJ) 和 最 大 简约 法 ( MP) 推断 这 些 种 的 系统 发 育 关系 。 通 过 Kishino-Hasegawa ( KH) 和 Shimodaira- 
Hasegawa ( SH) 检验 评估 这 4 种 系统 树 间 的 差异 ,确定 最 为 合理 的 系统 发 育 树 。【 结果 】 本 研究 新 
测序 获得 20 种 ( 亚 种 ) 的 COT 和 ITS2 序列 的 长 度 范围 分 别 为 625 ~685 bp 和 300 ~559 bp, CO T 
和 ITS2 序列 在 库 蚁 属 成 对 蚁 种 间 的 遗传 距离 分 别 是 0.002 ~0.198 和 0.006 ~ 1.807, 3A 23 
个 种 ( 亚 种 ) 的 COT +ITS2 序列 的 ILD 检验 结果 显示 ,数据 集 具 有 不 相 容 性 ,因此 CO TD. € ITS2 序 
列 不 适用 于 这 些 种 ( 亚 种 ) 的 系统 发 育 研 究 。 经 KH 和 SH 检验 显示 ,4 种 系统 发 育 树 中 ,基于 COT 
序列 构建 的 BI 树 最 为 合理 ,而 基于 ITS2 序列 构建 的 MP 树 最 为 合理 。 基 于 CO ID 核 葵 酸 序 列 所 构 
建 的 BI 树 显示 , 除 幼小 库 蚊 Cx. infantulus fe X& 2f PE3X. Cx. brevipalpis 外 ,各 亚 属 间 成 员 聚 类 结果 与 
传统 的 形态 学 归 类 结果 吻合 ;路 蚂 亚 属 Lutiza Je 6533. E Æ, Barraudius AR 92 N È 33: E 4, Culex 内 ; 梅 
zE H Maillotia ze 3$ È 3X E J£ Nexoculex R 7j — 30; ÈRA EH Culiciomyia, £L Z 3X. E A 
Eumelanomyia $e 3 f& 33: LÆ Lophoceramyia 均 显 示 为 单 系 。 基 于 ITS2 序列 所 构建 的 MP 树 显示 ,各 
亚 属 间 和 种 ( 亚 种 ) 间 关系 混乱 。 【结论 ] 重建 的 分 子 系统 发 育 树 表明 库 蚁 亚 属 不 是 一 个 单 系 群 。 
在 重建 库 蚊 属 系统 发 育 关 系 时 , 相 较 于 ITS2 fe COI +ITS2,COT 是 更 为 理想 的 分 子 标记 。 本 研究 
构建 的 分 子 系统 发 育 关系 为 中 国库 归 属 中 各 亚 属 和 各 种 ( 亚 种 ) 之 间 的 亲缘 关系 的 研究 商定 了 基础 。 
关键 词 : BUR 库 蚁 属 ; 分 子 标记 ; COI; ITS2; RARR; 中 国 
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Abstract: [ Aim] To sequence the CO Į and ITS2 sequences of the Culex species (subspecies) with 
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distribution record in China, and to construct and discuss the molecular phylogenetic relationships among 
these species ( subspecies). [Methods] CO I and ITS2 sequences of a total of 20 Culex species 
(subspecies) were newly sequenced, and CO [ sequences of other 20 Culex species ( subspecies) and 
ITS2 sequences of other three Culex species ( subspecies) were downloaded from NCBI databases. The 
base composition, interspecific genetic distance and saturation analysis of CO [ of 40 species 
(subspecies) and ITS2 sequences of 23 species (subspecies) were analyzed, and the incongruence 
length difference (ILD) test was conducted for CO [ + ITS2 sequences. The phylogenetic relationships of 
these Culex species ( subspecies) were inferred using maximum likelihood ( ML), Bayesian interference 
(BI), neighbor joining ( NJ) and maximum parsimony (MP) based on the nucleotide sequences of two 
molecular data sets (CO I and ITS2) , respectively. Kishino-Hasegawa ( KH) and Shimodaira-Hasegawa 
(SH) tests were used to evaluate the differences of these four phylogenetic trees constructed, and the 
most reasonable phylogenetic tree was selected. [ Results] The newly sequenced CO [| and ITS2 
sequences are 625 — 685 and 300 -559 bp in length, respectively. The pairwise genetic distances among 
Culex mosquito species ( subspecies) range from 0.002 to 0. 198 for CO I and from 0. 006 to 1.807 for 
ITS2. The ILD test of CO I + ITS2 sequences of 23 Culex species (subspecies) showed that the data set 
was incompatible, so CO | + ITS2 sequences were inappropriate to construct the phylogenetic trees of 
these species ( subspecies). By KH and SH tests among these four phylogenetic trees, the BI tree based 
on CO [| sequence was the most reasonable tree, while the MP tree based on ITS2 sequence was the most 
reasonable tree. The BI tree constructed based on the CO [ nucleotide sequences showed that all 
members, except for Cx. infantulus and Cx. brevipalpis of the same subgenus, are clustered together, 
which is consistent with the traditional morphological classification. Both the subgenus Lutiza and 
Barraudius are grouped into the subgenus Culex. The subgenus Maillotia and subgenus Nexoculex are 
grouped together. The subgenus Culiciomyia , subgenus Eumelanomyia and subgenus Lophoceramyia are 
monophyletic. The MP tree based on the ITS2 sequences showed that the relationships between subgenera 
and species (subspecies) are chaotic. [ Conclusion] The molecular phylogenetic tree reconstructed shows 
that the subgenus Culex is not a monophyletic group. Compared with ITS2 and CO Į - ITS2, COT isa 
more desirable molecular marker for the reconstruction of the phylogeny of Culex species. The molecular 
phylogenetic relationships reconstructed in this study lay a foundation for studying the phylogenetic 
relationships among subgenera and species (subspecies) of the genus Culex in China. 


Key words: Culicidae; Culex; molecular marker; CO [ ; ITS2; phylogeny; China 











库 蚊 属 Culex X je T Eb H 4j ( Insecta ) 390388 H 有 20 


属 38 种 (Poinar et al., 2000) 。 目 前 虽 已 有 一 


(Diptera ) 蚊 科 (Culicidae ) 库 蚊 亚 科 (Culicinae ) JÆ IX 
族 (Culicini) 。 全 世界 已 知 库 蚊 属 蚊 虫 26 亚 属 768 
种 ( 亚 种 )( Harbach, 2011) ,它们 广 布 于 世界 所 有 








些 库 尽 种 类 之 间 的 系统 发 育 学 研究 (Miller et al., 
1996; Demari-Silva et al., 2011, 2015; Harbach et 
al., 2012; Luo et al., 2016; Sun et al., 2017) ,但 对 











动物 区 系 。 中 国 已 知 库 蚊 属 蚊虫 9 亚 属 84 fh OO 
PR) (MEER, 1997; 董 学 书 等 , 2010, 2014) ,在 全 
国 范围 内 均 有 分 布 ( 付 文 博 和 陈 斌 , 2015 ) 。 库 蚊 属 
是 库 蚊 亚 科 中 的 第 一 大 属 ( Harbaech , 2011) ,它们 
是 一 类 重要 的 医学 昆虫 ,除了 吸血 和 骚扰 外 ,其 中 的 
一 些 种 类 还 是 班 氏 丝 虫 病 和 流行 性 乙 型 脑 炎 等 的 主 
要 传播 媒介 ( 陆 宝 鹿 和 赵 彤 言 , 2000) 。 蚊 科 的 基础 
分 类 研究 已 历时 200 余年 ,形成 了 较为 成 熟 的 体系 
( 师 永 霞 等 , 2008) 。 目 前 发 现 并 描述 的 蚊虫 化 石 仅 
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整个 库 蚊 属 的 系统 发 育 关系 仍 知 之 甚 少 。 

线粒体 DNA (mitochondrial DNA, mtDNA) 是 动 
物体 内 的 核 外 遗传 信息 载体 , 因 其 分 子 量 小 ,拷贝 数 
高 ,基因 组 中 一 般 无 内 含 子 、 转 座 因 子 和 间隔 序列 ， 
几乎 不 发 后 倒 位 . 易 位 等 畸变 与 重组 等 特点 , 常 被 用 
于 分 子 进化 研究 ( 陈 星 等 , 2012) 。 在 线粒体 基因 组 
的 37 个 编码 基因 (2 个 rRNA 基因 、13 个 蛋白质 编 
码 基因 和 22 个 tRNA 基因 ) 和 一 个 控制 区 中 ,由 于 
CO I 基因 进化 速率 适中 ,能 保证 足够 变异 的 同时 又 
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能 被 通用 引物 扩 增 ,因此 它 是 目前 用 于 物种 鉴定 的 
条 形 码 基 因 , 也 被 广泛 用 于 昆虫 系统 发 育 研 究 ( 王 
备 新 和 杨 莲 芳 , 2002) 。 

核糖 体 DNA (ribosome DNA, rDNA) 是 核 基因 , 
具有 多 拷贝 ,由 18S$, 5. 8S 和 28S RNA 编码 区 .基因 
间隔 (ICS) 第 1 内 转录 间隔 区 (ITS1 ) 和 第 2 内 转 
录 间 隔 区 (ITS2) 和 外 转录 间隔 区 ( EST) 组 成 。 在 前 
体 RNA 形成 ?rRNA 时 , ITSI 和 ITS2 会 被 剪断 ,不 参 
与 核糖 体 的 形成 ,因此 该 区 域 所 受 选择 压力 小 ,进化 
速度 快 , 常 被 用 于 近 缘 种 的 杀 缘 关系 研究 (Mukha 
and Wiegmann, 2002) ,ITS2 也 被 广泛 应 用 于 亲缘 种 
的 鉴别 (Collins and Paskewitz, 1996) 。 

对 于 不 同 阶 元 和 类 和 群 ,基于 适合 的 分 子 标记 可 
客观 地 反映 真实 的 亲缘 关系 。 近 年 的 研究 表明 , 作 
为 分 子 标记 的 CO T f ITS2 序列 具有 客观 的 进化 信 
息 特征 , 据 此 重建 的 系统 发 育 关系 能 较 大 程度 地 符 
合 真 实 的 种 间 关 系 ( Miller et al., 1996; REAS, 
2003; 师 永 霞 等 , 2008; Wang et al., 2012; Laurito et 
al., 2013; Chan et al., 2014; Khoshdel-Nezamiha et 
al., 2016) ,并 有 助 于 鉴定 形态 相似 的 蚊 种 。 

本 研究 基于 重庆 师范 大 学 昆虫 与 分 子 生 物 学 研 
究 所 的 馆藏 标本 ,鉴定 和 测序 了 库 蚊 属 20 个 种 ( 亚 
种 ) 的 CO I 和 ITS2 序列 ,同时 从 NCBI 数据 库 
(https: // www. ncbi. nlm. nih. gov/) 下 载 了 早 前 报道 
的 在 中 国 具 有 分 布 记 录 的 库 蚊 种 ( 亚 种 ) 的 CO E fil 
ITS2 序列 ,共计 获得 库 蚊 属 40 个 种 ( 亚 种 ) 的 CO T 



































( 亚 种 ) 都 具有 COT 基因 序列 ,其 中 的 23 个 种 ( 亚 
种 ) 具 有 ITS2 序列 。 其 中 ,20 个 种 ( 亚 种 ) 的 COI 
和 ITS2 序列 由 本 研究 测定 ,其 余 的 COT 和 ITS2 JF 
列 下 载 自 NCBI 数据 库 , 各 序列 来 源 见 表 1。 本 研究 
所 用 样本 于 2013 - 2016 年 采集 自我 国 广西 .山东 、 
黑龙 江 和 内 蒙古 ( 表 1) ,JE-T Bii BE (1997 ) 和 董 学 
书 等 (2010) 作 形态 学 鉴定 并 做 好 采集 记录 后 , 置 于 
-20% 冰箱 中 保存 ,以 备 CO I 和 了 ITS2 测序 。 

1.2 目的 片段 的 扩 增 及 测序 

按照 闫 振 天 等 (2017 ) 的 方法 ,每 种 蚊虫 选取 其 
1 头 蚊 虫 的 单条 后 足 提 取 基 因 组 DNA ,用 ND-1000 
全 波长 微量 测 浓度 仪 检测 其 纯度 和 浓度 后 ,保存 于 
-80% 箱 内 备用 。 使 用 Bio-Rad T100™PCR 扩 增 仪 
作 CO 和 了 ITS2 片段 的 PCR 扩 增 ,反应 体系 为 50 
uL, H} DNA 模板 1~2 pL, 正 反 引物 (10 pmol/L) 
各 1 uL, 2 x Taq MasterMix 25 nkL, 加 入 灭 菌 水 补足 
至 50 uL, 

用 于 CO L 基因 扩 增 的 引物 (Folmer et al., 
1994) 为 :LCO1490 (forward) :5'-GGTCAACAAATCA 
TAAAGATATTGG-3'; HC02198 (reverse) :5'-TAAAC 
TICAGGGTGACCAAAAAATCA-3'。 反 应 条 件 : 95*C 
预 变性 5 min; 95% 变性 1 min, 55?C 退火 1 min, 
72% 延伸 90 s, 35 个 循环 ; 72% 延伸 7 mins HF 
ITS2 扩 增 的 引物 (Beebe and Saul, 1995) 为 : 正 向 引 
物 ITS2A; 5'-TGTGAACTGCAGGACACAT-3'; Jz [5] 
引物 ITS2B; 5'-TATGCTTAAATTCAGGGGGT-3', Jx 

















序列 ,其 中 23 个 种 ( 亚 种 ) 具 有 ITS2 序列 。 本 研究 
进一步 分 析 了 这 些 序列 的 基本 特征 ,基于 这 些 CO I 
和 ITS2 序列 ,使 用 最 大 似 然 法 (maximum likelihood , 
ML)、 贝 叶 斯 法 (Bayesian inference，BI)、 令 接 法 
(neighbor joining, NJ ) 和 最 大 简约 法 ( maximum 
parsimony, MP) ) 分别 推断 中 国库 蚊 属 部 分 蚊虫 种 间 
的 系统 发 育 关系 , 并 使 用 PAUP * 4. 0 bl0 
( Swofford , 2002) 对 各 种 方法 构建 的 系统 发 育 树 进 
行 了 KH (Kishino-Hasegawa ) 检 验 和 SH (Shimodaira- 
Hasegawa) 检验 ,以 选择 出 最 合理 的 发 育 树 来 讨论 中 
国库 蚊 属 部 分 蚊虫 种 间 的 系统 发 育 关系 ,这 对 中 
库 蚊 属 昆虫 的 正确 鉴定 和 系统 发 育 关系 的 认 知 具 有 
科学 意义 。 
































1 材料 与 方法 


1.1 样本 和 序列 
本 研究 包括 了 库 蚊 属 40 个 种 ( 亚 种 ) ,所 有 种 








应 条 件 : 95% 预 变性 5 min; 95% 变性 40 s, 55'C iB 
火 40 s, 72% 延伸 1 min, 35 个 循环 ; 72% 延伸 6 

将 PCR 产物 上 样 至 1% 的 琼脂 糖 凝 胶 中 进行 电 
泳 , 然 后 用 EB 染色 ,在 紫外 灯 下 观察 并 切 下 目的 条 
带 , 用 DNA 胶 回 收 试剂 盒 ( Qiagen QIAquick PCR 
Purification Kit) 对 目的 条 带 进行 纯化 后 , 送 往 成 都 
垫 科 梓 申 生 物 公 司 用 ABI3730 基因 分 析 仪 进行 双 
向 测序 。 
1.3 序列 分 析 

先 将 双向 测序 结果 通过 DNAMANx 软件 拼接 ， 
同时 根据 测序 峰值 进行 人 工 校对 。 再 将 导出 的 一 致 
序列 输入 NCBI 网 站 ,进行 Blast 相似 性 检索 ,显示 
20 种 ( 亚 种 ) 库 蚊 属 蚊虫 的 CO 工序 列 分 别 与 基因 库 
中 对 应 的 种 ( 亚 种 ) 的 CO 工序 列 的 同 源 性 为 97% ~ 
100% ;20 种 ( 亚 种 ) 库 蚊 属 蚊虫 的 ITS2 序列 分 别 与 
基因 库 中 现 有 库 蚊 属 蚊虫 的 ITS2 序列 的 同 源 性 为 
97% ~ 99% , 验证 了 形态 鉴定 的 准确 性 。 使 用 
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MEGA 6.0 ( Tamura et al., 2013) 软件 ,对 3 组 分 子 
数据 集 (CO I , ITS2 和 CO I - ITS2 序列 ) 的 碱 基 组 


成 ( nucleotide composition )、 变 异 位 点 (variable 








sites) 简约 信息 位 点 (parsimonious sites ) 等 进行 分 
Tr ,并 以 Kimura-2a (K2P) 为 模型 ,计算 种 间 成 对 遗 
传 距离 。 









































































































































































































































表 1 本 研究 所 用 样本 和 序列 信息 
Table1 Information of samples and sequences in this study 
采集 时 间 序列 号 参考 文献 
亚 属 种 / 亚 种 采集 地 (年 /月 /日 ) GenBank accession no. References 
Subgenus Species/ subspecies Collecting site Collecting date 
(year/month/day) COI ITS2 COI ITS2 
Tp 区 小 库 收 蒙古 本 研究 本 研究 
EUER PIENE por 2013/8/6 MF278807 — MF288851 iir iin 
Barraudius Culex modestus Inner Mongolia This study This study 
zs peii Shaikevich «mui 
库 蚊 亚 属 Culex THEE KM233149 。 KU056508 pid. RU 
Culex pipien pipiens et al., 2016 Unpublished 
3e fA Pas PS 本 研究 本 研究 
RCM A 2015/7/14  MF278809 — MF288837 aud A 
Culex pipiens pallens Shandong This study This study 
卷 库 收 研究 CELA 
siia 2015/7/30  MF278808 ^ MF288853 zm died 
Culex quinquefasciatus Guangxi This study This study 
2 yk E E 出 月 Miller et al., 
RR KJ012240 U33039 TE TEAS 
Culex torrentium Unpublished 1996 
zie 8 AE Y 本 研究 本 研究 
EAM 2015/7/30 —— MF278812 — MF288848 Maa je 
Culex pseudovishnui Guangxi This study This study 
环 带 库 收 本 研究 本 研究 
TUAE 2015/8/1 MF278814 — MF288842 sk D 
Culex annulus Guangxi This study This study 
=H s pd 本 研究 本 研究 
2015/8/19  MF278813 — MF288846 io itd 
Culex tritaeniorhynchus Guangxi This study This study 
^1 78 p dy 本 研究 本 研究 
mici 2015/8/19  MF278815 — MF288836 adus ma 
Culex gelidus Guangxi This study This study 
As E Pg i 本 研究 本 研究 
TERR 2015/8/10 MF278816 — MF288849 id aid 
Culex theileri Guangxi This study This study 
— HE ge pd 本 研究 本 研究 
PARF 2016/7/17 MF278817  MF288839 i ma 
Culex bitaeniorhynchus Guangxi This study This study 
e a E hy 本 研究 ARS 
odd: 2015/8/19  MF278818 — MF288854 Hm uisi 
Culex fuscocephala Guangxi This study This study 
以 态 库 收 本 研究 本 研究 
WEFR | 2015/7/30 MF278819  MF288845 a "BR 
Culex mimeticus Guangxi This study This study 
INIA E 本 研究 本 研究 
下 2015/7/25 MF278820 MF288850 UB d 
Culex mimulus Guangxi This study This study 
KMA E 本 研究 本 研究 
2015/7/30 — MF278821 — MF288843 ig diu 
Culex murrelli Guangxi This study This study 
ERT Pts EL 本 研究 本 研究 
M i 2013/7/30 MF278822 — MF288838 xd 1s 
Culex orientalis Heilongjiang This study This study 
^E P hy Ki t al., 
fa B x DQ149239 umar et a E 
Culex hutchinsoni 2007 
迷走 库 蚊 TE Maekawa et al., 
Culex vagans 2016 
Magda Kinare al; 
id dl DQ310144 : UE 2 
Culex sitiens 2007 
^R Pe thy Ki t al., 
AARE DQ154167 - icol z 
Culex whitmorei 2007 
2 rs UR Hg M i 





Culex infula 


Unpublished 
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续 表 1 Table 1 continued 






















































































































































































采集 时 间 序列 号 参考 文献 
亚 属 种 / 亚 种 采集 地 (年 /月 /日 ) GenBank accession no. References 
Subgenus Species/ subspecies Collecting site Collecting date 
year niontb/day) COI ITS2 COI ITS2 
库 状 蚊 亚 黑 点 库 收 :研究 研究 
Pm x 2015/7/30 MF278811 — MF288841 Qe ue 
Culiciomyia Culex nigropunctatus Guangxi This study This study 
^1 gf] JA by x 研 究 研究 
HERE 2015/7/25  MF278810  MF288835 Qu ne 
Culex pallidothorax Guangxi This study This study 
琉球 库 疏 i 
MERE AB738139 i 未 出 版 " 
Culex ryukyensis Unpublished 
EE E ihi 
a LC104335 - TM - 
Culex kyotoensis Unpublished 
HANIF n Bp 6 本 研究 -研究 
ROLE Nc i 2016/1/8 MF278826 — MP288852 本 研究 TER 
Eumelanomyia Culex foliatus Guangxi This study This study 
短 须 库 蚊 D0424961 N Kumar et al., E 
Culex brevipalpis 2007 
b se pe thy K l. 
SHE DQ149238 = E < 
Culex malayi 2007 
pd He db «mu 
中 绳 库 蚊 a 未 出 版 
Culex okinawae Unpublished 
Ey i] -研究 研究 
T 2016/7/21 MF278825 — MF288844 Ad nid 
Culex hayashii Guangxi This study This study 
Ax > AE 355 ped < 出 
T t BOUE 带 纹 库 蚊 pan 未 出 版 
Lophoceraomyia Culex cinctellus Unpublished 
£T Jf] JA vy K t al., 
DN AY729981 z ed - 
Culex rubithoracis 2007 
小 库 收 Ki t al., 
Aiit DQ267691 x N E 
Culex infantulus 2007 
ZINE fg y K t al., 
小 型 库 蚊 EU259304 " umar et a 
Culex minor 2007 
须 E iiy zt H 
ARE AB738260 - T * 
Culex bicornutus Unpublished 
BAF ; m -研究 “研究 
DEPAR A 2015/7/30 MF278823 — MF288840 aT dd 
Lutzia Lutzia fuscanus Guangxi This study This study 
贪 y “研究 “研究 
ARMEN 2015/7/30 MF278824 — MF288847 本 研究 本 研究 
Lutzia halifaxia Guangxi This study This study 
ior V. EE ibi JAE nb zm 
iplc dett pue i 未 出 版 
Maillotia Culex hortensis Unpublished 
E, RE iey F H 58 pz iy < 出 月 
新 库 蚊 亚 必 边 留 库 蚊 oem — d 未 出 版 i 
Neoculex Culex trubensis Unpublished 
惊 骇 库 蚁 «m Miller et al., 
PAREN KR688733 U33036 TUA M 
Culex territans Unpublished 1996 
1.4 系统 发 育 树 的 构建 ILD (incongruence length difference) 检验 ,检验 联合 
本 研究 进行 系统 发 育 分 析 时 分 别 依据 COT, ”分 析 的 数据 是 否 具 有 相 容 性 。 
ITS2 和 CO I «ITS2 的 核 苷 酸 序 列 。 理 论 上 , 单 基 在 构建 系统 发 育 树 之 前 ,应 对 用 于 构建 系统 发 





因 对 于 系统 发 育 关 系 重建 的 信息 贡献 量 小 于 联合 数 。”” 育 树 的 分 子 数据 进行 碱 基 饱 和 度 分 析 , 以 判断 其 是 
d GME, 2015) ,但 由 于 不 同 基因 具有 不 同 的 进 ” 否 适用 于 系统 发 育 分 析 。 本 研究 利用 MEGA 6.0 $k 
化 速率 ,因此 在 对 联合 数据 分 析 前 ,需要 对 其 进行 MF, K2P 模型 分 别 计算 分 子 数 据 集 的 转换 值 
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(Ts) AHAH ( Tv) 和 K2P 距离 (Ts + Tv) ,以 转换 值 
+ Bit (EL SP Ah , TP 5 (E A E [EC RI] 4) 8 2 2 
图 ,以 评估 序列 的 碱 基 蔡 换 饱 和 度 。 若 图 形 呈 明显 
的 线性 关系 , 则 说 明 未 饱和 ,可 用 于 系统 发 育 分 析 。 

分 别 使 用 PAUP x 4. 0 b10, MrBays3. 1. 2 
( Huelsenbeek and Ronquist F, 2001) 和 MEGA 6.0 
软件 重建 各 蚊 种 之 间 的 系统 发 育 关 系 。 在 整个 分 析 
过 程 中 ,都 选用 冈比亚 按 蚊 An. gambiae (CO 工序 列 
的 登录 号 为 : DQ465316. 1; ITS2 序列 的 登录 号 为 : 
AF470229. 1) 作为 外 群 。 

4r 31H] NJ 法 和 MP 法 在 MEGA 6.0 软件 包 中 
构建 NJ 树 和 MP 树 。 在 构建 ML 树 之 前 ,利用 
Modeltest 3.7 ( Posada and Crandal, 1998) , PAUP x 
4.0 b10 和 MrMTgui 进行 联合 分 析 , 得 到 hIRT 和 
AIC (akaike information criterion ) 标准 选 出 的 最 优 模 
型 和 此 模型 基于 分 析 数 据 的 最 适 参数 ,二 者 得 出 的 
标准 不 同时 , 优先 选择 AIC 标准 下 的 模型 (Posada 
and Buckley, 2004 )。 以 AIC 为 最 佳 标准 ,选择 出 
CO T fll CO I +ITS2 的 核 昔 酸 序 列 的 最 优 模型 为 
GTR +G +1 模型 ;ITS2 的 核 背 酸 序列 的 最 优 模型 为 
HKY +G 模型 。 在 PAUP* 4.0 b10 中 ,对 不 同 的 序 
列 选择 对 应 的 最 优 模型 并 赋予 最 适 参 数 ,启发 式 搜 
索 ML 树 , 由 于 计算 量 较 大 ,对 ML 树 进行 500 次 自 
举 检 验 ; 在 MrBays3.1.2 中 ,采用 贝 叶 斯 推断 法 , 模 
型 选择 和 ML 法 相同 ,运行 100 万 代 , 每 隔 100 代 就 
取样 一 次 并 保存 ,最 后 舍弃 前 2 500 代 的 老化 树 ,对 
剩余 树 合 为 一 树 , 即 得 到 BI 树 。 对 上 述 4 种 方法 构 
建 的 系统 发 育 树 在 PAUP * 4.0 b10 中 进行 KH 检验 
和 SH 检验 ,比较 各 系统 发 育 树 间 的 差异 是 否 显 车 ， 
对 差异 不 显著 的 系统 发 育 树 构建 50% 合意 树 ( 吴 静 
等 , 2010), 
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2.1 CO1 的 序列 特征 

结果 显示 ,20 种 ( 亚 种 ) 库 蚊 属 蚊虫 的 CO 工序 
列 长 度 在 625 ~ 685 bp 之 间 。 通 过 和 下 载 自 NCBI 
数据 库 中 的 CO 工序 列 校对 排 齐 后 ,得 到 40 条 长 度 
均 为 523 bp 的 序列 ,无 终止 密码 子 的 出 现 , 无 碱 基 
的 插入 和 缺失。 其 中 变异 位 点 160 个 , 占 总 位 点 数 
的 30. 696 ;简约 信息 位 点 102 个 , 占 总 位 点 数 的 
19.5% 。 碱 基 A, T, G 和 C 的 总 平均 含量 分 别 为 
28.9% , 40.296 , 14. 896 和 16.2% ,A + 平均 含量 
为 69.1% ,明显 高 于 G - C 的 平均 含量 (30.9% ) 。 





















































其 中 密码 子 不 同位 点 上 的 碱 基 也 存在 明显 的 偏向 
性 ,A +T 平 均 含量 在 第 1 位 点 为 69.0% ,第 2 位 点 
为 71.9% ,第 3 位 点 为 66.2% 。 

40 种 ( 亚 种 ) 库 蚊 属 蚊虫 的 CO 工序 列 的 种 间 遗 
传 距离 变化 范围 在 0.002 ~ 0. 198 之 间 。 种 间 的 最 
大 遗传 距离 为 三 带 只 库 蚊 Cx. tritaeniorhynchus 和 短 
须 库 蚊 Cx. brevipalpis 之 间 ,达到 了 0. 198 ;种 间 最 小 
遗传 距离 0. 002 出 现在 尖 音 库 蚊 Cx. p. pipiens 和 
SUB Fg Y. Cx. quinquefasciatus 之 间 。 
2.2 ITS2 的 序列 特征 

20 种 ( 亚 种 ) 库 蚊 属 蚊虫 的 ITS2 序列 长 度 在 
300 ~ 559 bp 之 间 。 和 下 载 自 NCBI 数据 库 中 的 
ITS2 序列 校对 排 齐 后 显示 ,23 条 ITS2 的 核 苷 酸 序 
列 对 比 位 点 共 276 bp ,存在 大 量 的 碱 基 的 插入 和 缺 
失 。 其 中 变异 位 点 220 个 , 占 总 位 点 数 的 79.7% ; 
简约 信息 位 点 186 个 , 占 总 位 点 数 的 67.4% 。 碱 基 
A, T, G 和 fc 的 总 平均 含量 分 别 为 22.7%, 
23.896, 25.3% 和 28.296, A + 了 平均 含量 为 
46. 5% , G +C 的 平均 含量 为 53.5% 。 

23 种 ( 亚 种 ) 库 蚊 属 蚊虫 的 ITS2 序列 的 种 间 遗 
传 上 距离 变化 范围 在 0. 006 ~1. 807 之 间 。 种 间 的 最 
大 遗传 距离 为 1. 087, 均 为 叶片 库 蚊 Cx. foliatus 与 
贪 食 路 蚊 Lt. fuscanus AERE PES. Cx. fuscocephala 
之 间 ; 种 间 最 小 遗传 距离 (0.006) 则 出 现在 伪 杂 钱 
库 蚊 Cx. pseudovishnui PI $Æ Jy JÆ IX Cx. orientalis 
之 间 。 
2. 3 COI +ITS2 的 序列 特征 

对 获得 的 CO I AN ITS2 数据 进行 组 合 , 连 接 后 
得 到 CO I + ITS2 联合 基因 该 组 的 核 苷 酸 比 对 位 点 
共 740 个 ,其 中 变异 位 点 330 个 , 占 总 位 点 数 的 
44. 696 ;简约 信息 位 点 268 个 , 占 总 位 点 数 的 
38.6% 。 碱 基 A, T, G 和 C 的 总 平均 含量 分 别 为 
27.396, 35.596 , 18.2% 和 19.0% ,A + 了 T 平 均 含 量 
为 62.8% ,G +C 的 平均 含量 为 37.2% 。 对 23 个 种 
(ERR) CO I + ITS2 联合 数据 进行 ILD 检验 ,结果 
显示 数据 具有 不 相 容 性 (P=0.001), 因 此 CO T + 
ITS2 联合 数据 不 适用 于 进行 这 些 种 ( 亚 种 ) 的 系统 
发 育 研 究 。 
2.4 COI 和 ITS2 序列 数据 集 的 碱 基 饱 和 度 分 析 

2 组 分 子 数据 集 的 碱 基 饱 和 度 见 图 1。 图 1 
(A) 和 (B) 中 得 到 的 Ts 和 Tv 值 与 K2P 距离 之 间 均 
呈 明 显 的 线性 关系 ,这 说 明 在 遗传 进化 过 程 中 ,2 组 
分 子 数据 集 (CO 工序 列 和 ITS2 序列 ) 的 碱 基 的 转换 
和 其 换 均 未 达到 饱和 ,可 以 用 于 系统 发 育 分 析 。 
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图 1 不 同 K2P HER F CO T (A) 和 ITS2(B) 序 列 的 转换 (Ts) 和 颠 换 (Tv) 值 
Fig.1 Values of transitions (Ts) and transversions (Tv) for CO I (A) and ITS2 (B) sequences plotted against K2P distance 


2.5 库 蚊 属 40 TR (ERR) 的 系统 发 育 关系 

2.5.4 基于 COI 序 列 的 系统 发 育 关系 构建 :通过 
SH 和 KH 检验 对 基于 40 种 ( 亚 种 ) 库 蚊 属 蚊虫 的 
CO 工序 列 所 构建 的 4 种 系统 发 育 树 进行 了 评 佑 , 结 
RER BI 树 似 然 值 最 小 , 且 与 ML 树 、NJ 树 和 MP 
树 似 然 值 差异 显著 (已 <0. 05 ) , 即 相 较 于 其 他 3 种 
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此 ,本 研究 在 基于 CO 工序 列 所 构建 的 4 种 系统 发 育 
树 来 讨论 中 国库 蚊 属 蚊虫 的 系统 发 育 时 ,以 BI 树 为 
主要 依据 。 图 2 是 用 40 种 ( 亚 种 ) 库 蚊 属 蚊虫 的 
CO 工序 列 所 构建 的 BI 树 。 在 BI 树 中 ,各 亚 属 内 的 
种 ( 亚 种 ) 首 移 聚 类 , 且 与 形态 学 归 类 结果 ( 陆 宝 麟 ， 
1997) 基本 吻合 ,符合 库 蚊 亚 属 Culex BOVE Lutiza , 








构 树 方法 ,BI 法 所 构建 的 系统 发 育 树 最 为 合理 。 因 
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京都 库 蚊 Cx. kyotoensis LC104335 
琉球 库 蚁 Cx. ryukyensis AB738139 


42 99 mE Cx. bicornutus AB738260 
— E 小 型 库 蚁 Cx. minor EU259304 


73 带 纹 库 蚊 Cx. cinctellus AB738191 
红 胸 库 蚊 Cx. rubithoracis AY729981 


67 T5 叶片 库 蚁 Cx. foliates MF278826 
60 林 氏 库 蚊 Cx. hayashii MF278825 
马 来 库 蚊 Cx. malayi DQ149238 


库 状 蚊 亚 属 Culiciomyia、 新 库 蚊 亚 属 Nexoculex 、 绝 大 部 


尖 音 库 蚊 Cx. p. pipiens KM233194 
SÉ PESE Cx. quinquefasciatus MF278808 “| 库 蚊 亚 属 Culex 
HEER Cx. p. pallens MF278809 

凶 小 库 蚊 Cx. modestus MF278807 

类 迷走 库 蚊 Cx. torrentium KJ012240 

迷走 库 蚊 Cx. vagans LC054504 

希 氏 库 蚊 Cx. theileri MF278816 

AEE Cx. hutchinsoni DQ149239 

环 带 库 蚊 Cx. annulus MF278814 

WARREN Cx. pseudovishnui MF278812 
三 带 吃 库 蚊 Cx. tritaeniorhynchus MF278813 
白 霜 库 蚊 Cx. whitmorei DQ154167.1 

白雪 库 蚊 Cx. gelidus MF278815 

海滨 库 蚊 Cx. sitiens DQ310144.1 

棕 头 库 蚊 Cx. fuscocephala MF278818 

类 拟态 库 蚊 Cx. murrelli MF278821 

小 拟态 库 蚊 Cx. mimulus MF278820 

拟态 库 蚊 Cx. mimeticus MF278819 
KRAE Cx. orientalis MF278822 
Kmk Cx. infula AY729983 
ZE Cx. bitaeniorhynchus MF278817 
褐 尾 路 蚊 Lt. fuscanus MF278823 

SURE Lt. halifaxia MF278824 





| Gy CE JR Barraudius 


库 蚊 亚 属 Culex 





[psc Lutzia 


幼小 库 蚊 Cx. infantulus DQ267691 EAREE Lophoceramyia 
白 胸 库 蚊 Cx. pallidothorax MF278810 
黑 点 库 蚊 Cx. nigropunctatus MF278811 库 状 蚊 亚 属 Culiciomyia 











簇 角 蚊 亚 属 Lophoceramyia 











真 黑 蚊 亚 属 Eumelanomyia 
冲绳 库 蚊 Cx. okinawae AB738154 
98 惊 骇 库 蚊 Cx. territans KR688733 亚 属 Neocul 
边 留 库 蚊 Cx. rubensis LC054489 [roe P Neocilex 
霍 顿 库 蚁 Cx. hortensis KJ012080 | 梅 蚊 亚 属 Maillotia 
短 须 库 蚊 Cx. brevipalpis DQ424961 | 真 黑 蚊 亚 属 Eumelanomyia 
冈比亚 按 蚊 An. gambiae DQ465316 | 外 群 Outgroup 











图 2 基于 COT 核 芭 酸 序列 构建 的 40 种 ( 亚 种 ) 库 蚊 属 蚊虫 的 BEI 


Fig. 2 BI tree of 40 Culex species (subspecies) based on mitochondrial CO Į nucleotide sequences 
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分 真 黑 蚊 亚 属 Eumelanomyia HU $k fA tx MV JE 
Lophoceramyia 的 聚 类 , 且 置 信 值 均 大 于 5096 ; 库 状 
BOLJE 真 黑 蚊 亚 属 和 簇 角 蚊 亚 属 均 能 聚 为 单独 的 

支 , 置 信 值 均 大 于 60% ; 包 蚊 亚 属 Barraudius 和 路 
蚊 亚 属 归 人 库 蚊 亚 属 内 ,置信 值 都 大 于 80% ; 梅 蚊 
亚 属 Maillotia 和 新 库 蚊 亚 属 聚 为 一 文 , 置信 值 是 
55% 。 但 幼小 库 蚊 Cx. infantulus 和 短 须 库 蚊 在 BI 
树 中 的 分 类 地 位 与 传统 的 形态 学 分 类 结果 ( 陆 宝 



































ITS2 序列 所 构建 的 4 种 系统 发 育 树 进行 了 评估 , 结 
果 显 示 MP 树 似 然 值 最 小 , 且 与 BI 树 、ML BERI NJ 
树 似 然 值 差异 显著 (已 <0. 05 ) , 即 相 较 于 其 他 3 种 
构 树 方法 ,MP 法 所 构建 的 系统 发 育 树 最 为 合理 。 
因此 ,本 研究 在 基于 ITS2 序列 所 构建 的 4 种 系统 发 
育 树 来 讨论 中 国库 蚊 属 蚊虫 的 系统 发 育 时 ,以 MP 
树 为 主要 依据 。 用 23 种 ( 亚 种 ) 库 蚊 属 蚊虫 的 ITS2 
序列 所 构建 的 MP 树 见 图 3,MP 树 中 仅 符 合 库 收 亚 









































Wi, 1997) 差异 较 大 。 在 BI 树 中 ,属于 簇 角 蚊 亚 属 
内 幼小 库 蚊 与 库 状 蚊 亚 属 这 一 文 和 库 蚊 亚 属 ` 包 蚊 
亚 属 和 路 蚊 亚 属 组 成 的 一 大 支 互 为 并 系 关 系 ; 属 于 
真 黑 蚊 亚 属 内 的 短 须 库 蚊 则 是 和 男 外 两 支 ( 梅 蚊 亚 



































BRP RE E H (R E JE, Beha JE R E E 
Cx. p. pallens 和 类 迷走 库 蚊 Cx. torrentium ) , Y x. 
亚 属 和 库 状 蚊 亚 属 的 聚 类 ,置信 和 度 都 大 于 80% 。 与 
形态 学 分 类 绪 果 ( 陆 宝 刨 ，1997 ) 相 比 ,分 子 树 中 亚 























属 + 新 库 蚊 亚 属 这 2 个 亚 属 形成 的 一 支 和 库 蚊 亚 属 
+ 包 蚊 亚 属 + 路 蚊 亚 属 + 库 状 蚊 亚 属 + 真 黑 蚊 亚 属 
+ 艇 角 蚊 亚 属 这 6 个 亚 属 形成 的 一 大 支 ) 互 为 并 系 
2.5.2 基于 ITS2 序列 的 系统 发 育 关 系 构 建 :通过 
SH 和 KH 检验 对 基于 23 种 ( 亚 种 ) 库 蚊 属 蚊虫 的 






































届 间 和 种 间 关 系 较为 混乱 ,如 属于 真 黑 蚊 亚 属 中 的 
林 氏 库 蚊 Cx. hayashii 归 入 了 库 蚊 亚 属 内 ,同属 于 
库 蚊 亚 属 拟态 库 蚊 组 的 拟态 库 蚊 Cx，mimeticus、 东 
方 库 蚊 Cx. orientalis .小 拟态 库 蚊 Cx. mimulus 和 类 
WSE Cx. murrelli RA 3828 8| — 3c , HAT RUR: 
他 同 亚 属 或 不 同 亚 属 的 种 聚 类 。 
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淡色 库 蚊 Cx. p. pallens MF288837 
致 倦 库 蚊 Cx. quinquefasciatus MF288853 
类 迷走 库 蚊 Cx. torrentium U33039 

凶 小 库 蚊 Cx. modestus MF288851 

希 氏 库 蚊 Cx. theileri MF288849 

伪 杂 鳞 库 蚊 Cx. pseudovishnui MF278812 
类 拟态 库 蚊 Cx. murrelli MF278821 

东方 库 蚊 Cx. orientalis MF278822 
REEN Cx. hayashii MF278825 

环 带 库 蚊 Cx. annulus MF278814 
ZWE Cx. bitaeniorhynchus MF278817 
白雪 库 蚁 Cx. gelidus MF278815 

小 拟态 库 蚊 Cx. mimulus MF278820 
=E Cx. tritaeniorhynchus MF278813 
棕 头 库 蚊 Cx. fuscocephala MF278818 

白 胸 库 蚊 Cx. pallidothorax MF278810 

黑 点 库 蚊 Cx. nigropunctatus MF278811 
叶片 库 蚊 Cx. foliates MF278826 

惊 骇 库 蚊 Cx. territans KR688733 

拟态 库 蚊 Cx. mimeticus MF278819 
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图 3 基于 ITS2 序列 构建 的 23 种 ( 亚 种 ) 库 蚊 属 蚊虫 的 MP 树 
Fig. 3 MP tree of 23 Culex species (subspecies) based on ITS2 sequences 





CO 工序 列 的 K2P 距离 , tE HEAT TOKOE EFE SES 





3.1 库 蚊 属 的 种 间 遗 传 距离 


Hebert 等 (2003 ) 通过 统计 200 Fh eE H E E 

















目的 种 间 K2P 距离 是 大 于 0. 03 ,并 认为 这 种 方法 可 
以 推广 用 于 其 他 物种 的 种 的 鉴定 。Wang 等 (2012 ) 
采用 同样 的 思路 ,统计 了 中 国 15 属 122 种 蚊虫 
CO 工序 列 的 K2P 距离 ,结果 表明 ,超过 98% 的 蚊虫 
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遵循 种 间 K2P 距离 大 于 0. 02 ,种 内 K2P 距离 小 于 
0.02 的 规律 。 在 本 研究 中 , 除 尖 音 库 蚊 复 合体 ( 尖 
音 库 尽 、 致 倦 库 蚊 和 淡色 库 蚊 ) 内 的 种 ,其 他 库 习 种 
间 的 CO 工序 列 的 K2P 距离 均 大 于 0. 02 ,这 与 国内 
外 学 者 对 蚊 科 CO 工 基因 的 研究 结果 一 致 (Cywinska 
et al., 2006; Wang et al., 2012; Laurito et al., 
2013), 

在 本 研究 中 ,COI 基 因 的 K2P 距离 在 尖 音 库 蚊 、 
致 倦 库 蚊 和 淡色 库 蚊 之 间 的 变化 范围 为 0.002 ~ 
0.005。 相 似 的 情况 还 出 现在 另外 一 些 研究 中 , 例 
如 ,Laurito 等 (2013 ) 对 来 自 阿 根 廷 和 巴西 的 库 蚊 中 
的 致 倦 库 蚁 和 人 尖 音 库 蚊 的 遗传 距离 进行 分 析 发 现 ， 
二 者 之 间 的 遗传 距离 为 0.016 ,同样 小 于 0.02。 

关于 * 尖 首 复 合体 "是 一 个 多 态 性 物种 还 是 
妹 种 的 复合 体 一 直 具 有 争议 (Harbach，2012 ) 。 也 
有 一 些 支 持 后 者 观点 的 报道 ,比如 Liu 等 (2013 ) fl 
用 Ace-2 基因 研究 尖 音 库 收 复合 体 的 遗传 多 态 性 中 
就 证 明了 人 尖 音 库 蚊 、 致 倦 库 蚊 和 淡色 库 蚊 在 同一 区 
域内 可 以 相互 杂交 ;Cui S (2007) WEK T A F 
在 淡色 库 蚊 和 致 倦 库 蚊 的 杂交 体 。 因 此 ,有 关 尖 音 
库 蚊 、 致 伴 库 蚁 和 淡色 库 蚊 的 分 类 地 位 更 确切 的 结 
论 还 需要 进一步 调查 和 更 多 的 实验 及 数据 才能 
确定 。 
3.2 库 蚊 属 亚 属 和 种 的 系统 发 育 关 系 

本 研究 选择 了 2 个 分 子 标志 以 阐明 中 国 部 分 库 
蚊 属 蚊虫 的 亲缘 关系 ,分 别 为 CO I 和 ITS2 的 部 分 
序列 。 依 据 单 基 因 分 析 结 果 显 示 :使 用 ITS2 所 构建 
的 系统 发 育 树 中 ,除了 库 蚊 亚 属 中 尖 音 库 疏 组 ( 尖 
音 库 尽 、 致 倦 库 蚊 淡色 库 蚊 和 类 迷走 库 蚊 ) rg c 
亚 属 和 库 状 蚊 亚 属 内 的 各 成 员 能 聚 类 到 一 起 ,其 余 
的 亚 属 之 间 和 蚊 种 之 间 的 关系 都 较为 混乱 , 且 分 文 
的 置信 和 度 也 较 低 (小 于 50% ) 。 虽 然 ITS2 作为 分 子 
标志 在 解析 按 蚊 属 种 间 关 系 , 以 及 对 其 近 缘 种 和 隐 
藏 种 的 鉴定 时 具有 和 较 大 的 贡献 ( 马 雅 军 等 , 2011; 
马蜂 等 ,2013; Sum et al., 2014) ;但 是 将 其 运用 到 
库 蚊 属 蚊虫 中 ,可 能 不 能 得 到 较为 理想 的 结果 , 比如 
尖 音 库 蚊 复合 组 的 ITS2 序列 在 复合 组 内 亚 种 间 和 
亚 种 内 的 变异 存在 较 大 的 重 三 现象 ( 宋 社 吾 等 ， 
2003), 。 由 于 使 用 CO I 序列 所 构建 的 系统 发 育 关 
系 ,可 较 大 程度 地 客观 阐明 本 研究 库 蚊 种 间 亲 缘 关 
系 , 且 置信 值 大 多 高 于 50% 。 基 于 此 , 可 认为 在 本 
次 研究 的 库 蚊 属 部 分 蚊 种 中 , 相 较 于 ITS2 序列 ， 
CO 工序 列 更 适合 作为 重建 中 国库 蚊 属 蚊虫 系统 发 
育 关系 的 分 子 标 记 。 


































































































基于 形态 学 特征 所 建立 的 分 类 系统 中 , 相 较 于 
Kl; ac BE (1997 ) 的 分 类 系统 ,在 Harbach (2016) 的 新 
分 类 系统 里 ,传统 路 蚊 亚 属 Lutzia 升 为 库 蚊 族 内 路 
蚊 属 Lutzia ( Tanaka, 2003 ); 重 新 把 麻 蚊 亚 属 
Oculeomyia 从 传统 库 蚊 亚 属 内 的 二 带 唆 库 蚊 组 中 分 
离 出 来 (Tanaka, 2016)。 本 研究 中 ,基于 CO I 基因 
构建 的 BI 树 显 示 ,路 蚊 亚 属 和 二 带 嗓 库 蚊 组 中 的 二 
Tr UA Pg SY Cx. bitaeniorhynchus 都 聚 于 库 蚊 亚 属 内 ， 
且 置 信 值 都 大 于 50%., Wang % (2012) , Chan 等 
(2014) 和 Laurito 等 (2013 ) 同样 运用 CO 工 基 因 所 构 
建 的 系统 发 育 树 也 支持 路 蚊 亚 属 和 二 带 叭 库 蚊 组 聚 
于 库 蚊 亚 属 内 。 但 以 上 分 子 系统 发 生 树 均 依据 
CO [这 一 单 基因 所 构建 ,不 足以 涵盖 这 些 种 的 所 有 
遗传 信息 , 故 还 需 进 一 步 的 研究 来 精确 阐述 其 分 类 
地 位 。 重 建 的 系统 发 育 树 还 显示 , 库 蚊 亚 属 不 是 一 
个 单 系 群 ,这 一 结果 与 Harbach 等 (2012 ) 的 研究 结 
果 相 一 致 。 

本 研究 基于 CO I 和 ITS2 的 核 苷 酸 序列 ,使 用 
ML 法 .BI 法 \NJ 法 和 MP 法 重建 了 中 国 部 分 库 蚊 属 
蚊虫 的 系统 发 育 关 系 。 根 据 对 比 2 种 分 子 数 据 所 构 
建 的 系统 发 育 树 推测 :基于 CO I 基因 构建 的 BI 树 
更 适合 于 分 析 本 研究 中 所 涉及 的 库 蚊 属 蚊 虫 的 系统 
发 育 关系 , 且 CO 工 基因 是 分 析 库 蚊 属 系统 发 育 关系 
的 理想 分 子 标记 。 由 于 物种 有 限 和 所 用 分 子 标志 信 
息 量 的 限制 等 原因 ,并 不 能 精确 解析 库 蚊 属 内 部 各 
个 阶 元 的 系统 发 育 关系 。 但 本 研究 测定 的 20 种 
( 亚 种 ) 库 蚊 属 蚊虫 的 COT 和 1ITS2 序列 为 物种 鉴定 
和 中 国库 蚊 属 蚊 虫 的 进一步 系统 发 育 研 究 黄 定 了 一 
定 基 础 。 
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